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ABSTRACT 



Specification of parameters defining the mean, as well as the 
turbulent properties , of the over water surface layer is examined . A 
baroclinic boundary layer model is considered using Cardone's modifi- 
cations for the marine environment. 

Data fields of pressure and temperature were selected as input, 
and original ship observations were utilized to verify the predicted 
surface parameters. The effects of stability changes, isobaric curva- 
ture, thermal wind, and the value of the von Karman constant in the 
model were examined with results shown. 

Iteration criteria used in the model were evaluated with respect 
to operational application. Winds calculated by the model at the height 
of 19.5 meters were incorporated into the current FNWC wave height 
program, and the results compared with observed wave data. Turbulent 
parameters describing the properties of the index of refraction were 

also computed and evaluated with respect to laser propagation criteria. 



2 



TABLE OF CONTENTS 



I. INTRODUCTION 13 

n. ASPECTS OF A PLANETARY BOUNDARY LAYER 15 

A. SURFACE BOUNDARY LAYER 15 

1. Basic Structure of the Boundary Layer 15 

a. Viscous Sublayer 16 

b. Surface Layer 16 

c. Ekman Layer 17 

2 . Similarity Theory for the Surface 

Layer 17 

3. Consideration of the Drag Coefficient 

(C^) and Roughness Parameter (zq) 20 

B. PLANETARY BOUNDARY LAYER 24 

1 . Background 24 

2. Cardone's Extension of Blackadar's 

Two- Layer Model 25 

3 . Some Characteristics of the Marine 

Surface Boundary Layer 31 

III. COMPUTATIONS AND NUMERICAL PROCEDURES 35 

A. EXTERNAL WIND FIELD SPECIFICATIONS 35 

1. Geostrophic Wind Computations 35 

2. Correction for Curvature 37 

3. Thermal Wind Correction 38 



3 






■i 



J 






j 



i, 



B. PROCEDURES IN SOLVING THE BOUNDARY 

LAYER EQUATIONS 39 

1 . Initial Conditions 39 

a. Geostrophic Wind Speed Less Than 

or Equal to 5 . 0 meters/second 39 

b. Geostrophic Wind Speed Greater 
Than 5.0 meters/second and Absolute 

"' Value of Air-Sea Temperature Difference 
Less Than or Equal to 1°C 40 

c. Geostrophic Wind Speed Greater Than 

5.0 meters/second and Absolute 
Value of Air-Sea Temperature Difference 
Greater Than 1“C 40 

2 . The Modified Blackadar Baroclinic 

Boundary Layer Model 4 1 

C. COMPUTATION OF INDEX OF REFRACTION 

CONSTANT 43 

D. WAVE HEIGHT ANALYSIS 46 

IV. SPECIFICATION OF EXTERNAL PARAMETERS 48 

A. FNWC ANALYSIS 50 

B. SUBJECTIVE ANALYSIS 55 

V. RESULTS FROM EXPERIMENTS WITH THE MODEL 59 

A. COMPARISONS IN TERMS OF SURFACE ROSSBY 

NUMBER AND STABILITY CATEGORIES 59 

B. EXPERIMENTS WITH CHANGING THE PHYSICS 
IN THE TWO-LAYER MODEL AND CHANGING 
THE SPECIFICATION OF THE SURFACE 

GEOSTROPHIC WIND 64 

C. COMPARISONS OF COMPUTED WINDS WITH 

OBSERVED WINDS 78 



4 



I 



■■ 



1 ?': 



I 









D. RESULTS FROM WAVE HEIGHT COMPUTATIONS 82 

E. RESULTS FROM INDEX OF REFRACTION 

COMPUTATIONS 86 

VI. CONCLUSIONS 89 

APPENDIX A COMPUTER PROGRAM FOR MODEL 91 

BIBLIOGRAPHY 99 

INITIAL DISTRIBUTION LIST 102 

FORM DD 1473 104 



5 







1 



a 



■ i- ■ 



L: 







■■■ , I 



'i 



i 



• M 




3:'fK- 











i 



LIST OF FIGURES 



Figure Page 

1. Diagram of the Planetary Boundary Layer 16 

2. Various Suggested Forms of the Drag Coefficient 

(C 2 ) Versus Windspeed (U^) at 10 meters 23 

3. Wind Vector Diagram Applicable at Level h 26 

4. Variation of the Geostrophic Drag Coefficient with 
Surface Rossby Number and Dimensionless Stability 

Length (L*) (Cardone, 1969) 28 

5. Variation of Surface Cross Isobar Angle with Surface 

Rossby Number and Stability (Cardone, 1969) 28 

6. Geometrical Relationships Resulting from the 
Requirement of the Continuity of Wind Shear at the 

Internal Boundary Level h (Blackadar, 1965b) 29 

7. Reported Wind-Geo strophic Wind Ratio as a 
Function of Air-Sea Temperature Difference 

(Carstensen, 1967) 32 

8. Two-Layer Model Predictions of the Dependence 

of the Ratio of the 19 . 5 meter Wind Speed to the 
Geostrophic Wind Speed Upon Air-Sea Temperature 
Difference and Thermal Advection (Cardone, 1969) 32 

9. Subset of the FNWC Hemispheric Grid Used in this 
Study. (Note that all the following charts fit this 

scale exactly (1;30 million projection).) 36 

10. Surface Pressure (solid) and Air-Sea Temperature 
(dotted) Fields for 00 GMT 30 Nov 69 based on 

FNWC Data 51 

11. Geostrophic Wind Speed Analysis Based on Surface 

Pressure (Figure 10) 52 

12. Gradient Wind Speed Analysis Based on Surface 

Pressure (Figure 10) 53 



6 




i 










I 




4 

i 

I 



« 

'i 

1 









13. Percent Difference Field for Geostrophic-Gradient 

Wind Fields (Figure 11-Figure 12) 54 

14. Surface Pressure (solid) and Air-Sea Temperature 
(dotted) Fields for 00 GMT 30 Nov 69 Based on 

Plotted Data 56 

15. Geostrophic Wind Speed Field Based on Surface 

Pressure (Figure 14) 57 

16. Gradient Wind Speed Field Based on Surface 

Pressure (Figure 14) 58 

17. Same as Figure 4 Except Dotted Lines are Results 

from this Study 60 

18. Same as Figure 5 Except Dotted Lines are Results 

from this Study 61 

19. Geostrophic Drag Coefficient Versus Air-Sea 
Temperature Difference for Range of Geostrophic 

Wind 63 

20. 19.5 Meter Wind Speed Analysis (m. sec.^) for 

Case 1 (N/NC/NT) based on FNWC Data 66 

21. 19.5 Meter Wind Speed Analysis (m. sec. ^ ) for 

Case 2 (NN/NC/NT) Based on FNWC Data 67 

22. 19.5 Meter Wind Difference Field (%) for Case 2- 
Case 1. Positive % Indicate Increased Wind 

Speed Due to Stability 68 

23. 19.5 Meter Wind Speed Analysis (m. sec. ^) for 

Case 3 (NN/NC/T) Based on FNWC Data 69 

24. 19.5 Meter Wind Difference Field (%) for Case 2- 
Case 3. Negative Values Indicate Increased 

Wind Speed Due to Thermal Wind 70 

25 . 19.5 Meter Wind Speed Analysis (m . sec . ^ ) for 

Case 4 (NN/C/T) Based on FNWC Data 72 

26. 19.5 Meter Wind Difference Field (%) for Case 3- 
Case 4. Negative Values Indicate Increased Wind 

Speed Due to Curvature 73 



7 







4 




1= 



■■■■ 








v; 



I 




r ; 



I 



8 

1 



.; . 4 



Jt 



V - -1 



¥ 



h\ 



\ 



i 



27. 19.5 Meter Wind Speed Analysis (m. sec.^) for Case 5 

(NN/NC/T-k=0.35) Based on FNWC Data 74 

28. 19.5 Meter Wind Difference Field (%) for Case 3- 
Case 5. Negative Values Indicate Increased Wind 

Speed Due to Varying Value of k 75 

29. 19.5 Meter Wind Direction Difference Field (degrees) 
for Case 3 -Case 2. Positive Values Indicate 

Increasing or Veering Wind Direction 77 

30. Observed Surface Wind Speed (m. sec. Analysis 

Based on Plotted Data and is Subjectively Analyzed 79 

31. 19.5 Meter Wind Speed Analysis (m.sec.^) with 
All Correction and Specification Terms Included, 

Based on FNWC Data 80 

32. 19.5 Meter Wind Speed Analysis (m. sec.^) with All 
Correction and Specification Terms Included, Based 

on Plotted Data 81 

33. Computed Wave Height Analysis (feet) Based on 

FNWC Data for 00 GMT 30 Nov 69 83 

34. Computed Wave Height Analysis (feet) Based on 

Plotted Data 84 

» 

35. Observed Surface Wave Height Analysis (feet) Based 

on Plotted Data and Subjectively Analyzed 85 

o 

36. Distribution of Temperature Structure Constant (C^). 

Plotted Values Should Be Multiplied by 10“^ 87 



8 



\ 

4 

I 

I 




1 



{ 



TABLE OF SYMBOLS 




f 

G 

g 

H 

h 

K 

Kc 

k 

L 



L' 



Reciprocal of the turbulent Prandtl number (K^,/K ) 

ri m 

Boston constant = 0.81 

Blackadar constant = 0,0003 

Refractive index structure constant 

Specific heat of air at constant pressure 

Structure constant for temperature fluctuations 

Nondimensional drag coefficient at level z 

Rate of dissipation of turbulent energy 

Angle between the geostrophic wind and the 
thermal wind 

Coriolis parameter 

Surface geostrophic wind speed 

Gravity = 9.8 meters/ second 

Upward turbulent heat flux 

Height of surface layer 

Streamwise component of wave number 

Curvature 

Turbulent transfer coefficient for heat 
Turbulent transfer coefficient for momentum 
von Karman constant 

Lettau - Monin - Obukhov stability length 
Modified stability length 



9 



i 



1 

] 

j 

'I 

' I 

i 



i 

I 

? 



1 




1 



5 






.. i 

I 

' i 

I 



I 



Dimensionless stability length 

Measure of the "smearing" of temperature 
inhomogeneities 

Kinematic viscosity 

Baroclinic model parameters 

Dimensionless baroclinic model parameter 

Surface pres sure (millibars) 

Divergence terms 

Nondimensional temperature gradient 

One-dimensional spectrum for temperature 
fluctuation 

Nondimensional wind shear 
Integrated profile stability parameter 
Surface geostrophic wind angle 
Radius of curvature 
Gas constant for dry air 
Richardson number 

Mean density of surface boundary layer 

Dimensionless thermal wind 

Surface air temperature (°C) 

Sea surface temperature (°C) 

Potential temperature (°K) 

Scaling temperature 

Air temperature measured at level of 
thermometer (°K) 







^ I 




. f ‘ 



i 



Potential temperature at the surface (®K) 
Sea surface temperature (®K) 

Surface stress 
Friction velocity 
Wind at level h 
Surface wind 
Wind at height z 

u, V components of geostrophic wind 
u, V components of thermal wind 
Gradient wind 

Ageostrophic wind component 
Wave height (feet) 

Height above sea level (meters) 

Surface roughness parameter (meters) 



ACKNOWLEDGEMENT 



This research was conducted under the advice of 
Professor Kenneth Davidson, Department of Meteorology, with consider- 
able assistance from Dr. Jack Kaitala of the Fleet Numerical Weather 
Central. -Professor R. E. Elsberry reviewed the final manuscript. 

The project required utilization of the Computer Facility of the 
Naval Postgraduate School and environmental data from the Fleet 
Numerical Weather Central in Monterey, National Weather Records 
Center in Asheville, N. C. , and from participating units in the PARKA II 
experiment conducted in the Pacific in 1969. 

I would like to express my sincere thanks to my advisors and 
members of the above facilities and units , without whose help the study 
could not have been completed. Especially, I want to thank my wife, 
Evelyn, for her patience and understanding during this study and also 
for her typing and proofreading of the manuscript. 



12 




i 



t 






] 

1 






I. INTRODUCTION 



Parameterization of the boundary layer over the sea surface has 
been receiving increased attention in recent years. Fortunately much 
of the theory developed for low level turbulence over a land surface is 
also -applicable over the oceans. But the air-water interface does pre- 
sent some unique problems. For instance a fluid surface, rather than a 
fixed one, results in mobile surface roughness elements. Because of 
the presence of a fluid surface, any parameterization requires a knowl- 
edge of the processes which transfer heat and momentum energy across 
the boundary and ultimately influence such variables as wave generation 
and index of refraction. 

The purpose of this study was to examine and evaluate several 
features of a recently formulated marine boundary layer model which was 
designed for operational use. The model is a two-layer baroclinic model 
for the marine boundary layer and the formulation examined was that 
developed by Cardone (1969). Cardone based his formulation on an 
earlier two-layer model due to Blackadar (1965a). Features of interest 
in the model are related to the physics of the planetary boundary layer 
and include the effect of atmospheric stability, baroclinicity , coriolis 
acceleration and changes in the surface roughness parameter on the 
pressure driven surface wind. 
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The scope of the study was to examine the model in the context of 



its operational usage. As such, it was a case study in which the data 
represented pressure and temperature observations from sequential 
12-hour periods when a large cyclonic pressure system dominated the 
Central North Pacific Ocean. This particular synoptic situation provided 
the necessary relations between the scalar fields to examine features 
associated with baroclinicity , stability and wind speed categories. 

Because these data represent scalar fields which had to be de- 
fined by some analysis scheme before being used as input to the two- 
layer model, the differences which arise due to procedures in the external 
field specifications are considered. The effect of isobaric curvature on 
the surface wind fields is examined and this effect is compared to the 
changes associated with the physics in the model. The differences 
which may arise solely because the fields are defined by an objective 
versus a subjective analysis are considered. 

Finally, expressions describing the boundary layer are, by their 
very nature, empirical. As such, they contain empirical constants and 
the von Karman constant (k) is one whose value has recently been ques- 
tioned. The importance of a suggested change from k=0.40 to 0.35 was 
considered . 
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II. ASPECTS OF A PIANETARY BOUNDARY LAYER 



A. SURFACE BOUNDARY LAYER 

1 , Basic Structure of the Boundary Laver 

Recent advances in the study of atmospheric turbulence 
have led to a greater confidence in describing the distribution of wind 
and temperature in the atmospheric boundary layer over a solid surface. 
This progress has arisen largely through the application of the Monin- 
Obukhov similarity theory which has been supported by recent observa- 
tional studies. Cardone's contribution was primarily to extend 
Blackadar's two-layer neutral baroclinic boundary layer model over fixed 
terrain, and then to apply this model to the boundary layer over a fluid 
surface. The major difference between these boundaries is that the sea 
surface changes its character as the wind blows over it. 

In general, Blackadar's model treats the boundary layer as 
being formed by both turbulent friction and the coriolis acceleration. All 
processes are considered to be three dimensional, and the layer is char- 
acterized by a density stratification. It is conveniently separated into 
the three regions depicted in Figure 1. In this delineation the important 
consideration is the specification of Kj^ which is the turbulent transfer 
coefficient for momentum. 



15 






( 

1 









< 600 m 



t 

z 




Free Atmosphere 



Ekman Layer 



Surface Layer 
^Viscous Sublayer 






<50 m 



< 1 cm 



K, 



m 



Fig. 1. Diagram of the planetary boundary layer 



a. Viscous Sublayer (less than 1 cm) 

In this near surface layer the stress is supported by 
viscosity (aerodynamically smooth flow), and the mean velocity 
increases linearly according to 



U 

z 



U* z 



( 1 ) 



where v is the kinematic viscosity and is the friction velocity. 

b. Surface Layer (up to 50 meters) 

This is the constant stress layer in which the coriolis 

parameter can be neglected and K increases linearly with height. For 

m 

neutral stratification the logarithmic law is valid 
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( 2 ) 



- u* 

U = “TT" log z 
z k o 

where k=0.4 is the von Karma n constant. 

c. Ekman Layer (50-600 meters) 

This layer is characterized by the stress decreasing 

with height and K constant throughout the layer, 
m 

In formulating the boundary layer model, Cardone (1969) 
considered the validity of the "log law" for the surface layer over water, 
which has been experimentally documented down to near surface limits 
(Phillips, 1966; Paulson, 1967). He also assumed that the process of 
wave generation did not seriously affect the surface layer since under 
neutral and steady conditions a logarithmic shear zone should be estab- 
lished. Accepting the log profile as being valid, Cardone was able to 
apply much of the surface layer turbulence theory developed over land 
to over-water conditions . However Cardone did take into account the 
action of the waves when specifying the effective roughness parameter 
for the sea surface. 

2 . Similarity Theory for the Surface Layer 

The Monin-Obukhov "similarity theory" has led to consider- 
able progress in the specification of the distribution of wind and tem- 
perature in the boundary layer. The basis of this theory is that near 
the ground there exist velocity (U^) , length (L), and temperature (0^) 
scales which are essentially invariant with height. When the principle 
variables such as temperature (T) , wind (U), and height (z) are expressed 
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as nondimensional fractions of these quantities, a series of nondimen- 
sional equations can be formed which are of general validity in the sur- 
face boundary layer (Lumley and Panofsky, 1964). Recent analyses of 
the oceanic wind data (e.g. Phelps and Pond, 1971) have indicated that 
the Monin-Obukhov theory is valid for over-water application. These 
results support Cardone's inclusion of the Monin-Obukhov stability 
criteria in the boundary layer model. 

Similarity theory predicts that a universal relation should 
exist for the nondimensional wind shear 



M _ w 

U* d z u 



(3) 



and the nondimensional temperature gradient 



?L. ^ ^ 

9* 3z 



^t 



(4) 



where and 0^ are functions of the dimensionless ratio z/L. Here 
0 is the potential temperature, 9* is the scaling temperature defined by 



9 = 1 lU- 

* kU* Cp Pa 

and L is the Lettau-Monin- Obukhov stability length defined by 



(5) 



. _ -U^ Cp PaT ... 

^ ~ kgH • 

The heat flux (H) is not easily measured, so a modified stability length 
was defined by Panofsky (1963) 

U*T hxx/bz 



L' 



kg 39/dz 



(7) 
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The assumption that L’ is independent of height implies that a, is con- 

h 

stant, which in turn implies similarity between wind and temperature 
profiles - 

If L* is a valid scaling length, similarity theory predicts 



that 



= 0 (z/L') 
0* ^z aj^ 



Integration of these formulae (Panofsky, 1963) yields 

U* 






0=0 + 

z o ah 



l^ln z/Zq - V(z/L')j 
- j^ln z/Zq - ip(z/L')j , 



where the relation between and 0 ^ is 



4>(z/L') = f 

n 



2/L' 1-^^ ii) 

I 



d^ , 1 = z/L'. 



( 8 ) 



(9) 



( 10 ) 



( 11 ) 



( 12 ) 



Paulson (1970), using the KEYPS formula, integrated (12) by parts to 
obtain the following expression for 



= 1 - - 3 In + 2 In (L^) + 2 tan" +ln (^).(13) 



1+ffi 



Within the scope of an operational application, these 
expressions for the surface layer would be used to determine the wind 
speed given U* , or to determine U* given the wind speed at a particular 
level. For this reason the partial derivatives in (7) are replaced by 
expressions involving height differences and air-sea temperature 
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differences. This is achieved by substituting the following expression 

for , obtained by combining (9) and (11) 

o z 

z 



ae _ ^u ~ Qs ) 

Sz ■ z [ln(^a/zo)- H>(^VL')] 



where z^ is the height at which is measured (chosen as 10 meters by 
Cardone) . 



Substituting (14) and (8) into (7) yields 

Z; 



L' = 



uf 9 [in (-^) - H^(^a^L')] 



k^g (0a - Qg) 

Equation (10) is then solved for U* 

U* = 



(15) 



k U 



m 



ln(^m/z^)_ v*;( 2 m/L-) 



(16) 



where z^ is the height at which the wind speed is measured (chosen as 
19.5 meters by Cardone). Equations (15) and (16) can then be solved 
for U* given the wind speed at some level (U 2 ) and the air-sea tempera- 
ture difference . 

3 . Consideration of the Drag Coefficient ( and Roughness 
Parameter (^o) 

The surface roughness parameter and the influence of thermal 
stratification have received considerable attention in the definition of 
the relation between the nondimensional drag coefficient (C 2 ) and the 
wind speed at a certain height (usually 10 meters). The need for such a 
relation arises from the "bulk aerodynamic" formula for estimating 
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momentum flux. This equation, deduced from dimensional analysis, has 



the form 



T = U2 



C U2 



(17) 



where T is the surface stress. From this expression the drag coefficient 
is defined as 



C 



z 




(18) 



z 

Rearranging and squaring (10) yields the following expression for the 
drag coefficient at a level z in terms of and stability 



t2 



C_ = 



In 



Vzq - y\>r/v) 



(19) 



which was the form considered by Cardone. 

Observations have indicated that over the sea, z^ is depend- 
ent on U*. This dependence, however, changes as the wind speed 
changes. For the case of light winds (aerodynamically smooth flow), 
the drag coefficient C_ decreases with increasing wind speed (Kraus, 
1966) and Zq has the form 



For aerodynamically rough flow, 
speed according to 



V 

1.09 U* ' 

increases with increasing 

tt2 

a U* 

g 



( 20 ) 

wind 



( 21 ) 



where (a) is a constant = 0.035 (Charnock, 1955). A summary of field 



* 










♦ 



observations indicates that the change between rough and smooth flow 
occurs at wind speeds of 5-6 meters/second. 

Cardone proposed the following simple expression for z 

o 

covering a wide range of wind speeds (both rough and smooth flow) 

^1 2 

Zq = — + C 2U* + C3 . (22) 

The constants were chosen so that was a minimum for a value of 
6 meters/second at the 10 meter height under neutral conditions, and 
above this speed the relation approximated Charnock's results. For U* 
in meters/second and z^ in meters, the resulting expression was 

Zq = 6.84 X 10 " Vu* + 4,28 X 10"^U^ - 4.43 x lO"."^ (23) 

Other representations of C versus U have been published. 
A few of these appear in Figure 2 along with Cardone' s approximation. 

It is clear from the scatter of curves in Figure 2 that the specification 
of the drag coefficient was a primary decision in Cardone' s formulation. 
This choice of C„ versus U_ will be examined during the discussion on 
results from this investigation. 

With the equations for stability length (15), friction velocity 
(16), and surface roughness parameter (23), the surface layer wind dis- 
tribution is completely described within the framework of the similarity 
theory. These three equations can be solved simultaneously for U* 
using iteration techniques with the external parameters described below 
as input. 
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Fig. 2. Various suggested foras of the drag 

coefficient (C^) versus wlndspeed (U^) 
at 10 meters 
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B. 



PLANETARY BOUNDARY LAYER 



1 . Background 

In developing a model for the entire boundary layer, con- 
sideration has to be given to the possible external parameters. In 
Cardone's formulation these parameters were the geostrophic wind (G), 
the air-sea temperature difference ), the coriolis parameter (f), 

and the thermal wind (r). He considered previously suggested empirical 
relations between the geostrophic drag coefficient (U*/G) and the sur- 
face Rossby number (Lettau, 1959) 



R = 
o f 



( 24 ) 



where G is the surface geostrophic wind speed. Empirical relations 
have also been suggested between the cross isobar flow angle and 

^o* 

Blackadar (1965a) developed a two-layer boundary layer 

\ 

model, consisting of a Prandtl surface layer and an Ekman layer, which 
was verified on the basis of the above relations. 

Cardone's formulation was closely related to Blackadar' s 
work. In particular, the height of the surface layer (h) was specified 
explicitly in terms of the external parameters 

Bo G 



h = 



(25) 



.-4 



where Bq was a nondimensional constant equal to 3. Ox 10 as deter- 
mined by Blackadar. 



24 



a 



i 



♦ 



By applying proper boundary conditions and requiring that 
continuity of wind and wind shear exist at level h, the following solu- 
tions for U*/G and (along with the solution for z^) completely 

define the surface boundary layer wind distribution under neutral 
conditions . 



Zi. 

G 

U* 



3/2 



[2 k sin^ <4*„] 

kV? sin (n/4- 4>o) 



Ro 



(26) 



(27) 



Implicit in such an approach is that the value of K at the internal 

m 



boundary between the surface layer and the Ekman layer is given by 



K 

m 



kU* G 



f 



(28) 



2 . Cardone's Extension of Blackadar's Two-Layer Model 

The success of Blackadar’s two-layer model in predicting 
the essential characteristics of the neutral surface boundary layer, along 
with its relative ease of application, suggested its extension to the non- 
neutral marine boundary layer, Cardone extended the diabatic model to 
over- water flow by regarding Zq as a function of U* . 

The eddy viscosity in the constant stress layer for the non- 
neutral case was expressed by Cardone in terms of the nondimensional 
wind shear as 



kU*z 

■m = 

/u (VL-) 



(29) 
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Since Kj^ is constant with height above h, its value in the Ekman layer 



is then 



K = 



kU* BqG 






(30) 



The wind speed at h was given by 



U* 

U, = — 

h k 



l^ln - ip (VL-)j . 



(31) 



Internal boundary conditions applied at level h appear in Figure 3 and 



follow from the fact that 



always makes an angle of 



dz 



to 



the ageostrophic wind (W) and that 



d U 



is parallel to U in the 



constant stress layer, 




Fig. 3. V/ind vector diagram applicable at level h. 

By requiring continuity of wind shear at h, as well as directional 
continuity, Cardone was able to reduce the problem (analogous to (26) 
and (27)) to 

U* r 0 

— = 1^2 k Bq sin^ iVq 0 ^ (VL‘) J (32) 

U* = k-^ sin (n/4-4^o) 

G In BqRq - IV (h/L') 
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u* 

These equations are solved simultaneously for — and 4 ) from the 

G o 

external parameters G, f, and In addition, a dimensionless 



stability parameter was defined as 

L* = 



BqG 



f L' 



(34) 



which represented the ratio of the surface layer h to the stability 
length L' (Blackadar and Ching, 1965). 

The relationships determined by Cardone for U*/G and 
as functions of and L* appear in Figures 4 and 5. For all surface 
Rossby numbers , U*/G increases and decreases with decreasing 

stability. In general, stable conditions have a more pronounced effect 
on the surface layer wind characteristics than unstable conditions. In 
comparing the results from the model with observations , Cardone noted 
that U*/G increased by roughly 30% in unstable conditions, with no 
marked tendency to decrease at great instability, and decreased roughly 
70% in stable conditions. These features in the observed results will 
be discussed later. 

The effects of baroclinicity were also included by assuming 
that the geostrophic wind was a linear function of height (Blackadar, 
1965b). Blackadar (1965a) introduced a parameter, p, which was defined 
at level h by 



P 



^h 

OU/^z)j^ ■ i0^{h/V) 



In BqR^ -IP (VL‘) . 



(35) 
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Fig. 4. Variation of the geostrophlc drag coefficient 
with surface Rossby number and dimensionless 
stability length (I^) (Cardone, 1969) 




Fig. 5. Variation of surface cross Isobar angle with 
surface Rossby number and stability ^Cardone, 

1969) 
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The requirement of continuity of wind shear for the baroclinic case took 



the form 



[' 



U. - P K |f)H- 



aw 

a 2 



L 



(36) 



Figure 6 illustrates the geometrical implication of (36). From the Ekman 
layer solution of the ageostrophic wind, the following equation can be 
obtained 



aw 
a 2 



= p- I Wj^ I , 



(37) 



where p' is a dimensionless parameter defined by 



BoG 

' 0u(hA') 






0 ^ (Vl') 

2kU* B^G 
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(38) 



Fig, 6. Geometrical relationships resulting from the 
requirement of the continuity of wind shear at. 
the internal boundary level h (Blackadar, 1965b). 



Next the thermal wind vector was replaced by its dimensionless 



magnitude 

i_ I aG 

^ f I a 2 



(39) 
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where the direction of 



^ G 
d z 



relative to the surface geostrophic wind 



was implied by angle n Figure 6. Still another dimensionless 
quantity was defined 



r' = rB rinBR-m (VL‘)1 
o'- O O ' / /J 



(40) 



and (35) was then rewritten as 



hG 



3 z 



= r' G, 



(41) 



By applying the law of cosines to each of the two smaller triangles in 
Figure 6 that have M for one side, the following relations were determined 



= 


1 1 2 2 2 
|Wj^ 1 s = G q 


(42) 


2 

s = 


(1 2p' + 2p'^) 


(43) 


2 

q 


(1 r'2 - 2 r' cos p ). 


(44) 



n is the angle between the geostrophic wind and the thermal wind (see 
Figure 6). Then the law of sines was applied to yield 

sin i = P'/s (45) 

sin y = — sin H (46) 

q 



The following relations apply to the large triangle (ABC) in Figure 6, 

Uh = |Wj^ I sin ( a + 7 )/ sin 

2 






W, 



- 2G 



I cos ( a -I- 7 ) 



= uj + G^ - 2 U^G cos m 
n ho 



(47) 

(48) 

(49) 



From equations (31), (45), (46), (47), (48), and (49), the following 
equations describing the surface boundary layer characteristics were 



obtained 
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Tan ly 



sin ( a + 7 ) 



(50) 



o 



(s/q)- cos ( a + 7 ) 
k G q sin ( “ + 7 ) 



U 



* 



s sin W [In B^R - Ul (VL')] ' 



(51) 



By including (23) and (15), Cardone obtained a system of equations that 
could be solved for U* and iP ^ from the input parameters G, f, 



knowledge of sea level pressure, surface air temperature, sea surface 
temperature, and latitude. 

3 . Some Characteristics of the Marine Surface Boundary Laver 



well the planetary boundary layer model, using only large-scale synoptic 
parameters as input, specified the surface boundary layer wind distribu- 
tion over the sea surface. Cardone considered several studies relating 
the geostrophic wind ratio U^/G and the air-sea temperature difference. 
The results of an extensive study by Carstensen (1967) at Fleet Numerical 
Weather Central (FNWC) appear in Figure 7 along with earlier results 
due to Bleeker. The ratio of surface wind (V) to geostrophic wind (G) 
is tabulated for each degree of the reported air-sea temperature differ- 
ence, and the median value for each interval is plotted as an x in 
Figure 7. The validity of the results was limited to the air-sea tempera- 
ture difference range -4 to +1°C. In order to include observed winds 
in the surface pressure analysis, FNWC uses a constant value for V/G 
of 0.78 and a turning angle of 15° toward lower pressure. This selection 
could produce significant error. 




, and n . These equations were calculated from a 



A primary effort of Cardone 's work was to determine how 



31 



c 



♦ 



1 .00 



99 



0.90 - 



0.80 - 

V 

Vgs 

0.70 - 



0.60 - 



10 




s 



\ 




-7 -6 -5 -'4 -3 -’1 5 'i j ■■ 

Ta - Tg "C 



Fig. 7* Reported wind-geostrophlc v.'ind ratio as a 
function of air-sea temperature difference 
(Carstensen, 19^7) 
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Fig. 8. Two layer model predictions of lihe dependence 
of the ratio of the 19.5 meter wind speed to 
the geostrophlc wind speed upon air-sea 
temperature difference and thermal advection 
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Cardone's model can be used to compute a wind speed at 
19.5 meters. Figure 8 illustrates the variation of the ratio of the 19.5 
meter wind and surface geostrophic wind with the air-sea temperature 
difference. The ratio varied less in unstable conditions than in stable 
conditions . 



The effects of baroclinicity also appear in Figure 8 and 
these effects increase as stability decreases, with the ratio, U , 

increasing with cold advection and decreasing with warm advection. 
Mendenhall (1967) also studied the influence of baroclinicity and sta- 
bility on the veering of the wind in the planetary layer. He concluded 
that over mid-latitude oceans the diurnal variation of frictional veering 
of the wind with height was negligible, but that lapse rate and baro- 
clinicity were important. Results obtained by Cardone, using the plan- 
etary model, showed qualitative agreement with Mendenhall's conclu- 
sions. Cardone's results also suggested that baroclinicity was more 
significant than stability in determining the veering of the wind in the 
planetary boundary layer. For both cases in situations of strong, cold 
advection the actual wind may back with height, especially in conjunc- 
tion with unstable stratification. 

Additional results on the influence of baroclinicity came 

from a study by Clarke (1970), who found a tendency for the boundary 

layer wind to back instead of veer with height in unstable conditions . 

He stated that this could have been due to baroclinicity or to accelera- 

d T * 

tion. He noted that it could have occurred if ^ ^ was systematically 
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negative in convective conditions up to a level of about 0 . 2 — 

instead of being positive as in the normal case. This latter interpreta- 
tion was based on the observation that the wind appeared to be sub- 
geostrophic in the convective layer and reached geostrophic values only 
after a rapid increase through the overlying inversion. 
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III. COMPUTATIONS AND NUMERICAL PROCEDURES 



A. EXTERNAL WIND FIELD SPECIFICATIONS 

All computations were made with data defined on a subset of the 
FNWC 63X63 grid which covered the North Pacific Ocean (Figure 9). 

Eata fields included sea level pressure, sea surface temperature, and 
surface air temperature at each of 315 grid points spaced 381 km apart 
at 60 °N. The necessary map factor was included for a 1:30 million 
polar stereographic projection. 

The Coriolis term appears as an external parameter and therefore 
poses a limitation on the lowest geographical latitude where the geo- 
strophic relation would be valid. It has already been established that 
the dynamic boundary layer was formed by the surface stress and the 
coriolis acceleration. Pushistov (1970) concluded that manifestations 
of the nonlinearity in the dynamics of the boundary layer and the appear- 
ance of a perceptible ageostrophic component of velocity in the free 
atmosphere were to be expected near latitudes of 7-10°. Sheppard 
( 1970) suggested 20°N as a limiting latitude for application to boundary 
layer models. A limiting value of 15° latitude was used in this program 
as the limit on geostrophic flow. For latitudes below 15°, the coriolis 
parameter would be kept constant. 

1 . Geostrophic Wind Computations 

The geostrophic wind components were computed from the 
surface pressure field (P) as follows 



35 







) ^ 






I 




36 




I 









f 



i 
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'^g Pa f 

_i a p 
'^g P f ^ X 

using a centered explicit finite difference technique . 
determined from the equation of state 



(52) 

(53) 

The density was 



_ ^ P 

Pa R<j 

where is the gas constant for dry air, and 
ture in degrees Kelvin. The magnitude of the 
computed at each grid point from 



(54) 

0^ the surface air tempera- 
geostrophic wind was 




(55) 



2 . Correction for Curvature 

A gradient wind was calculated to correct for the effect of 
curvature on the surface geostrophic wind. Two basic principles were 
required for this correction. The first was that following an air parcel 
along an isobar, the total derivative was zero or 



dP(x,y) = 0. 



(56) 



The second requirement was the mathematical statement of curvature 

= [l + (y)2]3/2 ^r' 

where y’ and y" were the first and second derivatives with respect to x, 
and R was the radius of curvature. By expanding the total derivative of 
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P, obtaining y' and y" , and rearranging, the following result was ob- 
tained for the radius of curvature 



R = 



( 

^ y 


r 


1’" 


bp f 

b X 


. P b^ P aP/3x1 

^y by^ b 1/^y J 


^P ^^P 
1 ^y ^x2 



(58) 



Next the gradient wind was computed from 



Vg = G- 



V2 

G 

fR 



(59) 



Solving for V and rearranging led to a workable form of the equation 



= 



(60) 



G 1/2 + yj 1/4 + G/Rf 
For this case the magnitude of the term G/Rf determined how much curva- 
ture would be allowed. The maximum values were chosen to be -3/16 
for anticyclonic curvature and +3/4 for cyclonic curvature which yielded 
a range of = 1.33G for anticyclonic curvature and = 0.67G for 
cyclonic curvature. This agreed with generally accepted synoptic limits. 

Subroutine GRADWD (Appendix A) makes the necessary com- 
putations and follows a similar scheme utilized by FNWC (Kaitala, 1971 
unpublished). Centered finite difference schemes were used to calculate 
the derivatives . 

3 . Thermal Wind Correction 

The nondimensional thermal wind components were deter- 
mined from 

(61) 



^T 



f2 e 



^ y 
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(62) 



V, 



S 



f2 ^ hx ' 



and the total thermal wind was given by 



r 




(63) 



For this model, r had a limiting value of 100. The geostrophic and 
thermal wind angles were calculated from the previously computed wind 
components . This procedure was repeated for each grid point and was 
used as input data to the boundary layer model. 

B. PROCEDURES IN SOLVING THE BOUNDARY LAYER EQUATIONS 
1 . Initial Conditions 

a. Geostrophic Wind Speed Less Than Or Equal to 5.0 
meters/ second: 

This case was associated with the aerodynamically 
smooth flow region described earlier. For this case no corrections were 
made for stability and a constant drag coefficient of 0.022 was assumed. 
This resulted in a constant inflow angle of 15°. The 19.5 meter wind 
was determined from 



The choice of the drag coefficient was not arbitrary, but depended on the 

relationship between C and U (see Figure 2). The value of C^^ was 

z z Xu 

a minimum at 6 meters/second under neutral conditions. The value of 

C,_ - was extrapolated from this curve. 

X ^ • o 



u 



19.5 



0.7G 



(64) 
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b. Geostrophic Wind Speed Greater Than 5.0 meters/ 

second and Absolute Value of Air-Sea Temperature 

Difference Less Than or Equal to 1°C: 

Initially the drag coefficient was C = 0.0245, with 

z 

the thermal wind (r) and the angle between the geostrophic and thermal 
wind ( n ) set equal to zero. For this case L^, L' and ly (h/L') were 
set equal to zero and 0 ^ (h/L')= 1.0. A solution was then obtained using 
the modified baroclinic model of Blackadar described below. 

c. Geostrophic Wind Speed Greater Than 5.0 meters/ 
second and Absolute Value of Air-Sea Temperature 
Difference Greater than 1° C : 

The drag coefficient C^= 0.0245 is used to estimate 
the initial U* and r and r| are set equal to zero. 0- - 0 was restricted 

3 S 

to the range -15 to +4.0°C. The values of L^, L', Ul (h/L'), and 
0 (h/L') were computed from z , 0 -0 and U and a solution was 

M 3. 3. S 

obtained using the modified Blackadar model described below. 

L^ was computed from (34) and L' from (15) using as 
a first approximation 

iy(^a/L') = 10 - In (^a/zo) . (65) 

iy(h/L') was computed from SUBROUTINE PSI (Appendix A). If L^ was 
greater than zero (stable case), then 

iy(h/L') = 0.7 L^ (66) 



40 





I 



If was less than or equal to zero (unstable case), a value of 
Richardson number and (h/L’) had to be computed. This was accom- 
plished in SUBROUTINE SHR (Appendix A), Then U) (h/L') was calculated 
from (13). 

JZfy (h/L‘) was computed using SUBROUTINE SHR 
(Panofsky, 1968; Cardone, 1969). If was greater than zero (stable 
case), then 



jefy (h/L') = 1 + 7L^ (67) 

If L^ was less than or equal to zero (unstable case), then a new 
Richardson number was determined from 



where initially 
inequality 



Ri = L^ 
■^new * 

^iold " • 



1/4 

(l-18Rioid> (68) 

This iteration was continued until the 



Ri 

•‘■new 

was satisfied, where 
found from 




(69) 



e = 0.01 for this model. Then 0 ^ (h/L') was 



= TTini — r V4 • 

^new ^ 

2 . The Modified Blackadar Baroclinic Boundary Layer Model 
The basis of this model was an iteration procedure which 
included the effects of stability and baroclinicity on the computed 
wind. The procedure was free of computational instability and converged 
rapidly for all reasonable choices of initial U* . First values of Zq , 
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P, P'/S,q, a, 7, / and U were obtained using respec- 

new 

tively equations (23), (35), (38), (43), (44), (45), (46), (50), and (51). 
The iterative procedure was repeated until the inequality 



U* - U* < < 
new ~ 



(71) 



was satisfied where e was 0.05. The iteration was repeated, setting 
, until the inequality was again satisfied. If L. was not 






new 



zero, then new values of ly , , and were computed where the 

argument of was z^/L’ . If 



T ' - T ' 

^ new ^ 



> 5 



(72) 



then L' was set equal to . The model was repeated from where 

L* was first calculated until the inequality was satisfied. After com- 
pleting the iterations for the case where the thermal wind (r) and n 
were equal to zero, the computed values of r and r| were substituted 
and the model was again repeated until all criteria were satisfied. If 
L’?^0, thenJ^(z/L’) was computed as J^(19.5/L') ; otherwise 0 {19 ,5/L') 
was set equal to zero. Finally the wind at 19.5 meters was determined 
from 



U 



U. 



19.5 



In - UJ ( ^^ 



19. 5 x *1 

J 



(73) 



Values of 9 -9 , G, z , z/L’ , i*i , U , U ^ , log R , U /G, 

as o o^iy.o 

U-Q /G and Rj were then available for each grid point. 
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c. 



COMPUTATION OF INDEX OF REFRACTION CONSTANT 



In this section possible relationships are considered between the 
refractive index structure constant and quantities dependent only on the 
wind speed and air-sea temperature difference. Specifically those 
quantities that can be predicted from the model considered in this study 
are examined. 

The propagation of light waves and radio waves is affected pri- 
marily by air-density inhomogeneities and therefore by temperature fluc- 
tuations . One effect of these fluctuations is scintillation or rapid changes 
in intensity at a point. These are distortions caused by variations in the 
speed of propagation through turbulent elements having different densi- 
ties, and hence different indexes of refraction. The refractive index 
structure constant (0^^) is the primary parameter characterizing the turbu- 
lent fluctuations of the atmospheric refractive index for scales from a 

2 

few millimeters to several meters. If the quantity 0^^ is known, cal- 
culations can be made on scintillation, beam spreading, and beam wander, 
all of which are important in describing light wave propagation. 

A relation necessary in the formulation is the one-dimensional 
spectrum for temperature fluctuations in the inertial subrange described 
by Corrsin (1951) 




(K) = B C 



-Vs 




(75) 



where e is the rate of dissipation of turbulent kinetic energy, N is a 
measure of the "smearing" of temperature inhomogeneities, B is an 
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enpirical constant (=0.81, Boston, 1970), and K is the streamwise 

component of wave number. The structure constant for temperature fluc- 
2 

tuations (C^ ) is defined as 




-1 



B e 






N 



(76) 



and is related to the index of refraction as follows 



,-79x 10 ® 

^ 2 f 

n 



(77) 



a 

Another expression which has to be considered is the nondimen- 
sional balance equation for turbulent kinetic energy, which, for station- 
ary, horizontally homogeneous flow is (Busch and Panofsky, 1968) 



0 .. - Vl - - 0 ^ = 0 



u 



(78) 



0 ^ is the nondimensional wind shear (8), 0 -^ represents the diver- 



gence terms , and 



0r 



e k z 



U' 



(79) 



The following relations arise from the definitions in the non- 



dimensional turbulent kinetic energy balance and definition of N 

.3 



U 



C = 0 



N = -Q w 



C kz 
^0 



(80) 



. = K„ (-|-^) . 

o z H o z 



(81) 



Substitution of these two expressions, plus the definition of K (29) 



m 



into the definition of yields (Panofsky, 1968) 



Ct B 0 ^ 



-1/3 K 



H 



K 



m 



0a 



(kz) ( 3 ^) . 



(82) 
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Equation (82) is not a suitable expression because of the lack of 

information on . However if the effect of the divergence terms 

is neglected, 0 can be expressed in terms of 0 and z/L as 

u 



0^ = 0 - z/L . (83) 

E u 

It is noted that with wind-wave interaction, the assumption that the 

divergence terms can be neglected is probably not valid, but the possible 

2 

effect this assumption has cn predictions of has not been studied. 
Substitution of (83) into (82) yields 



C" = B 
T 
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^H 1 


1 

N 





38 2 

(kz) (-r — ) . (84) 

o z 



This final equation for lends itself to the boundary layer model be- 
cause all the parameters can be computed within the model. Thus, the 

2 

model can be used to obtain a horizontal distribution of over water 
for various wind and stability conditions. The average stability ^€/^z 
is determined from (14). 

The ratio of the exchange coefficients K„/K is set equal to 1.0 

n m 

in Cardone's formulation. However it has been shown that and 
are not usually identical except for small values of dO/ bz for which 
the turbulent motion is not significantly influenced by buoyancy forces . 
Busing er et al. (1971) suggested the following form for the ratio of eddy 
diffusions 

n _ Q z/j \ ^2 

(85) 



a = 1.35 , . 

h (1-15 z/L) V4 
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Here a^= 1.35 for neutral conditions, which agrees closely with labora- 
tory measurements (cf . , Hinze, 1959). A value of 1.35 rather than 1.0 

for a, under neutral conditions arises because the von Karman constant 
h 

was equal to 0.35 instead of the usual 0.40 (Businger et al. , 1971). 



D. WAVE HEIGHT ANALYSIS 

Wave analysis is another application for the boundary layer model. 
For this purpose, a simplified version of the FNWC sea-swell model 
(Hubert and Mendenhall, 1970) was utilized. The u and v components 
of the 19 . 5 meter wind determined from the boundary layer model were 
used as both current and history input in this analysis. 

The FNWC sea-swell model computes a weighted wind from the 
wind components as follows: 

u = 0.3 0-7 u^ 



V = 0.3v^_j^ 2 0*7v^ 



(86) 



where t-12 is the 12-hour old history wind speed and t is the current 
wind speed. These component values are saved as history for the next 
analysis run, and then the total wind (V) is computed. The wave heights 
(W^ ) in feet are then computed from 





<4V>' 




3.3 

) 



+ 1.0 



(87) 



Although this is a very simple empirical method for determining wave 
heights, the basic method has been used by FNWC since the mid- 1960' s 
and this particular model has been used operationally since August 1970. 
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In addition, the FNWC model computed corrections for fetch and 



ice, and performed some smoothing of the resulting height field. It was 
assumed these corrections would be minimal for this application. 
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IV. SPECIFICATION OF EXTERNAL PARAMETERS 



A primary objective of this study was to apply the two-level baro- 
clinic boundary layer model to actual atmospheric conditions and compare 
the computed results with the observed results. The time period 00 GMT 
28 November 69 through 00 GMT 02 December 69, with an analysis inter- 
val of every 12 hours, was chosen for this comparison for several reasons. 
First an intense low pressure system, which influenced most of the central 
North Pacific Ocean, developed during this period with the system reach- 
ing maximum intensity about 00 GMT 30 November 69 . The system pro- 
vided input data which allowed the model to predict boundary layer 
parameters under varying atmospheric conditions . Secondly considerable 
synoptic data had been gathered by the author during an earlier study on 
the extreme sea conditions (unpublished) associated with this storm. The 
latter gave a good basis for comparing computed values to observed data. 
Finally, this storm produced waves and swell which caused considerable 
damage to several islands of the Pacific as well as to the spar buoy FLIP. 
Therefore a wave generation model was applied using the predicted sur- 
face winds , and these predicted waves were compared to the observed 
wave conditions . 

The 00 GMT 30 November 69 period was of the most interest due 
to its extreme conditions, extensive data coverage, and it was the prob- 
able period when the destructive waves were being generated. The various 
corrections described above were incorporated into the model for this 
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time period, and synoptic data were plotted and analyzed for surface 
pressure, air-sea temperature difference, surface wind speed and wave 
height. 

Surface layer parameters were computed including various combina 
tions of stability, curvature, and thermal wind to yield information on 
the relative importance of each parameter in the formulation. The follow 
ing notation will be used in describing the results: 

NN - non-neutral, correction for stability 
N - neutral - no stability correction 
NC - no curvature correction applied 
C - curvature correction applied 

NT - no thermal wind correction applied 
T - thermal wind correction applied 

It was recognized at the onset of the study that the specification 
of the external fields ( A0 and G) had to be reasonably accurate in 
order to evaluate the model on the basis of observed wind fields. For 
example, if the observed surface winds were quite different than those 
predicted by the model, some estimate had to be available on the possi- 
bility that the discrepancies were due to poor specification of the geo- 
stophic wind or the air-sea temperature difference. Therefore in this 
section consideration is given to the relevant aspects of representing 
those fields which make up the external parameters . 

The charts in the following figures correspond to a 1:30 million 
polar stereographic projection with the grid as depicted in Figure 9. 
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A. FNWC ANALYSIS 

Sea level pressure and air-sea temperature difference fields for 
00 GMT 30 November appear in Figure 10. A 967 mb low pressure 
center was located near 45 °N 180° with considerable thermal instability 
both ahead of and behind the system. A significant feature of this sys- 
tem was the intense isobaric gradient in the western quadrant. 

The first consideration was the comparison of the computed surface 
geostrophic wind with the computed gradient wind for this period . The 
geostrophic and gradient wind fields computed from the surface pressure 
fields in Figure 10 appear in Figures 11 and 12 (units of meters/second). 
Figure 13 depicts the influence of the curvature correction on the geo- 
strophic wind computed as a percent difference from geostrophic wind 
minus gradient wind then divided by geostrophic wind. 

The two maximum wind centers west and north of the Low are based 
on data at one grid point and result from the influence of anticyclonic 
curvature, strong isobaric gradients, and finite difference computations. 
These features do not appear to be realistic. However, the wind speed 
distribution depicted by them is good except that the magnitude of the 
wind speed may be a little too large. It was expected that the maximum 
change would occur in an area of maximum curvature and maximum pres- 
sure gradient from the relation in (60). Generally speaking, the change 
in wind speed due to the curvature correction is 10 to 20% with a maxi- 
mum change of 34% imposed by the limits placed on the term G/Rf in (60). 
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B. SUBJECTIVE ANALYSIS 



The results from the hand analyzed surface chart were also used 
In comparisons of results. Figure 14 illustrates the surface pressure 
and air-sea temperature difference fields from the subjective analysis 
which are analogous to Figure 10. The major variations between Figure 
10 and 14 appear to be due to the difference in the number of synoptic 
reports used as well as the analysis scheme. Theoretically, Figure 14 
should be more representative since more data were included, but it was 
a subjective analysis. Figures 15 and 16 depict the geostrophic and 
gradient surface winds computed from the subjectively analyzed pressure 
field . 
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• Fig, l6. Gradient wind speed field based on surface 
pressure (Figure 1A) 
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V. RESULTS FROM EXPERIMENTS WITH THE MODEL 



A. COMPARISONS IN TERMS OF SURFACE ROSSBY NUMBER AND 
STABIUTY CATEGORIES 

Near surface parameters were computed for the 19.5 meter level. 

'This level represents an average height of ship anemometers (Bunting, 

1968). Theoretically the computed 19.5 meter level winds could be 

compared with the observed wind fields and this would yield some 

information on the accuracy of the model in predicting surface winds . 

In order to verify the computations and also to examine the range 

of conditions provided by the data, a comparison was made with Cardone's 

results which appear in Figures 4 and 5. Results for this comparison 

were obtained from computations in which stability was included but not 

baroclinicity . Results from all nine time periods appear in Figures 17 

and 18 where values of the logarithm of the surface Rossby number (24) 

range from 7.76 to 9.45. The U./G versus log R results (Figure 17) 

^ o 

agree with Cardone's results (Figure 4) except with regards to the inter- 
val between L^ isolines. In Figure 17 this interval (spacing) is larger. 
With respect to 4^^ versus log R^ (Figure 18), the unstable results 
(negative L*) agree with those in Figure 5 except, again, with regards 
to the interval between L* isolines . The stable results (positive L^ ) , 
however, deviate considerably from Cardone's results with the veering 
angle decreasing as L^ increases. The results in Figure 18 would 
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Fig. 17, Same as Figure 4 except dotted lines are results from this study 
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Fig, 18. Same as Figure 5 except dotted lines are results from this study 




t 



indicate that the veering angle decreases with both Increasing stability 
and increasing instability with a discontinuity occurring near neutral 
conditions . 

The unusual, and unexpected, features in the results just de- 
scribed were examined in considerable detail. The same results were 
noted by Pafias (unpublished study, 1971) in an investigation with 
Cardone's program. These features in the results are perhaps related 
to the drag coefficient representation selected by Cardone and also be- 
cause different iteration limits were used in this study than those used 
by Cardone . 

The range of data with respect to U./G and 0 -9 was also 

* as 

examined and the results appear in Figure 19. The computations for 
these results were made on the basis of a barotropic atmosphere with no 
correction for isobaric curvature. These computations differ from those 
used by Cardone for Figure 8 in that a variable coriolis parameter was 
used. Isolines have been drawn for geostrophic wind speeds ranging 
from 15 to 75 knots. For those wind speeds considered by Cardone in 
Figure 19, and also for wind speeds above 20 knots considered by 
Pafias, these results are in agreement. The change in shape of the 
curves for winds under 20 knots was not described by Cardone but is 
probably due to a combination of the functional forms for the stability 
dependence and the drag coefficient as chosen by him. 
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Fig. 19. Geostrophlc drag coefficient versus air-sea 

temperature difference for range of geostrophlc v;lnd 
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B. EXPERIMENTS WITH CHANGING THE PHYSICS IN THE TWO- LAYER 

MODEL AND CHANGING THE SPECIFICATION OF THE SURFACE 

GEOSTROPHIC WIND 

In this section ^ results which were obtained from different 

computations with the two-layer model are compared. These computations 
were varied by including or not including specific parameters or correc- 
tions . The variations considered and compared were described earlier 
and are 1) N/NC/NT, 2) NN/NC/NT, 3) NN/NC/T, 4) NN/C/T, and 
5) NN/NC/T with k=0.35. 

These variations will be referred to as Cases 1 through 5. It is 
noted that if the geostrophic wind speed was less than or equal to 5 . 0 
meters/second, stability and baroclinicity were not taken into account 
in the iteration schemes. In regions with this wind speed range, espe- 
cially in the subtropics, differences will not arise except where isobaric 
curvature exists . 

Per cent difference (1-2) will be used to compare the two cases 
where (1-2) corresponds to the following computation at each grid point: 

case 1 result - case 2 result . 

case 1 result ' 

Although large differences between cases were unexpected, small cor- 
rections at low wind speeds could produce large per cent changes. 

The surface wind fields corresponding to each of the five cases 
are presented in the following paragraphs along with selected comparisons 
between pairs of cases. 
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Case 1 (N/NC/NT) 

This is the simplest case with respect to the possible physics 
which are included in the model. The external wind field is defined for 
a neutral, barotropic atmosphere with no curvature correction. The wind 
field at 19.5 meters for this case appears in Figure 20. 

Ca se 2 (NN/NC/NT ) 

The effects of stability are now included and the 19 . 5 meter wind 
field for this case appears in Figure 21. The per cent difference (2-1) 
field is shown in Figure 22. Positive values correspond to increases in 
the computed U ^ field due to including the stability effect. This in- 
crease is quite evident on the western side of the cyclone. The patterns 
are not expected to be coincident with the air-water temperature isolines 
in Figure 10 because the stability effect is dependent on both the temper- 
ature gradient and the wind speed. The 20% contours near the eastern 
border are attributed to small corrections at low wind speeds as dis- 
cussed previously. 

Case 3 (NN/NC/T ) 

The effects of including the baroclinicity or thermal wind are now 
included and the 19.5 meter wind field for this case appears in Figure 
23. The per cent difference (2-3) field is seen in Figure 24. Negative 
values correspond to an increase in the computed 19.5 meter wind due 
to including the thermal wind in the model. In this difference field, only 
20 grid points (6%) appeared to be significantly affected by including the 
thermal wind. These grid points generally lie in an area of warm air 



advection. 
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5 meter wind speed analysis (m. sec 
e 1 (N/NC/NT) based on FNWC data 
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Fig. 21* 19*5 meter wind speed analysis (m. sec 
Case 2 (NN/NC/NT) based on FKWC data 
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Case 4 (NN/C/T) 

In addition to all specifications possible within the physics of the 
model, the correction for curvature is applied to the external wind field. 
The 19.5 meter wind analysis for this field is shown in Figure 25 and 
the per cent difference field (3-4) appears in Figure 26. Negative values 
correspond to larger 19.5 meter winds due to including the curvature cor- 
rection. Strong wind maximums, noted in comparing the external wind 
fields with and without this correction, also appear in Figure 26 as the 
predominant feature. In general, the effective change in the 19.5 meter 
wind field specification is larger and covers more area than any of the 
other cases . 

Case 5 (NN/NC/T with k=0.35 ) 

This case differs from case 3 only with respect to the change of 
the von Karman constant from 0.40 to 0.35. The 19.5 meter wind field 
for this case appears in Figure 27 and the per cent difference field (3-5) 
is depicted in Figure 28. Negative values correspond to increased 19.5 
meter winds due to the change in the value of k. The suggested value 
of k=0.35 arose from recent surface layer observational investigations 
by Businger et al. (1971). They noted from their results that using 
k=0.40 in several empirical expressions led to estimates of values 
which were 15% too high when compared with direct measurements of 
the momentum transfer. The fields in Figure 28 indicate that a general 
increase of 12 to 17% occurs in the computed 19.5 meter wind field if 
the lower value of k is used. This increase agrees with the 
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FIs. 28. 19.5 meter wind difference field (^) for Case 3 
Case 5. Negative values Indicate Increased wind 
speed due to varying value of k 
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discrepancies noted by Businger et al. . Because the von Karman 
constant appears in several expressions in the boundary layer model, 
the importance of its change in the two-layer model was not known. 
These results indicate that the difference between k= 0.35 and k=0.40 
is important. 

In summary, the five cases described above indicate which terms 
were most influential in the model. It should be remembered that most 
of the subtropical points were not influenced by the model since the 
wind speeds were generally less than 5.0 meters/second. Of all the 
modifications to the model, the curvature correction to the external 
field produced the largest variation over the entire grid. Changing k 
from 0.40 to 0.35 accounted for the second most significant change. 
Differences in the surface wind fields were evident, however, when 
specifications in the model were made. Of the two possible specifica- 
tions, stability appeared to be more significant than the thermal wind 
with respect to area influenced and per cent difference changes . 

Although the influence of the thermal wind in this model appears 
to be small with respect to the relative magnitude of the computed wind, 
the influence it has on the direction of the wind appears to be signifi- 
cant. The difference field, in degrees, of the wind direction for 
NN/NC/T - NN/NC/NT appears in Figure 29. This field represents the 
difference in wind direction between a case with and a case without the 
thermal wind specification. The region of significant differences does 
not coincide with that region where the thermal wind effect was large 
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with respect to wind speed. Positive values in Figure 29 correspond to 
a veering wind and negative values correspond to a backing wind. It is 
noted that the direction change due to baroclinicity is primarily deter- 
mined by the angle between the thermal and geostrophic winds . 

C. COMPARISONS OF COMPUTED WINDS WITH OBSERVED WINDS 

These comparisons will be described for the case NN/C/T with 
k=0.35. The 19.5 meter wind fields computed from the other cases were 
compared to the subjectively analyzed observed surface wind field 
(Figure 30) and in all cases the computed winds underestimate the ob- 
served winds. This was also the case when the above corrections were 
included, but yielded the best results with respect to the magnitude and 
patterns of the wind speeds. 

Because the previous comparisons revealed the importance of the 
specifications of the external wind field with respect to curvature, it 
was considered necessary in this comparison to examine the importance 
of the initial analysis (FNWC or subjective). The 19.5 meter wind field 
computed with the above specifications and with the FNWC analysis is 
shown in Figure 31 and the wind field based on subjective analysis is 
seen in Figure 32. The observed surface wind field which appears in 
Figure 30 was analyzed using only those observations reported above 
10 meters/second. In Figures 31 and 32, the maximum contour, 30 
meters/second, agrees with the observed maximum wind speed, but the 
objective analysis does not indicate the 25 meter/second maximum east 

of the low. 
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•Fig, 31 , 19,5 meter wind speed analysis (m, secT^ with all 
correction and specification terms included, based 
on FNWC data 
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A comparison of computed 19.5 meter winds from the two analyses 
(Figures 31 and 32) and the external wind fields (Figures 11 and 15) 
shows that differences in the external wind fields are consistently 
carried through the model. The model, however, was able to produce 
19.5 meter wind fields in both cases which were in good agreement with 
the observed wind field, at least for wind speeds over 10 meters/second. 

D. RESULTS FROM WAVE HEIGHT COMPUTATIONS 

Wave height fields were computed using the 19.5 meter wind 
fields obtained from the boundary layer model and the simplified wave 
height program developed from the FNWC operational sea-swell model. 
The program was run starting at 00 GMT 28 Nov in order to build history 
for this analysis. The 19.5 meter wind computations were obtained 
using all specifications and corrections described previously, and they 
were also made on both the FNWC and subjectively analyzed fields. The 
computed wave field from the FNWC data is seen in Figure 33 and the 
wave field from the subjective analysis appears in Figure 34. The ob- 
served wave field obtained from synoptic observations is depicted in 
Figure 35. It is noted that wave height is probably the most difficult 
parameter to observe and therefore a representative analysis of the ob- 
served wave field is difficult to achieve. 

The wave field obtained from the subjective analysis. Figure 34, 
appears to be in good agreement with the observed heights with maxi- 
mum computed heights slightly less than the maximum observed heights. 
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Results computed from the FNWC analysis, Figure 33, appear to have 



similar agreement with respect to the maximum height but poorer agree- 
ment with respect to the overall pattern of the wave field . 

In general these results indicate that the two- layer model can 
produce surface wind fields which can be used for wave height analysis, 
even from a very simple wave model. It is noted that in the subtropics 
the wave heights computed from both sets of analysis are considerably 
below the observed wave heights . This is probably due to low surface 
geostrophic wind speeds in this area and the fact that no corrections are 
made for stability or baroclinicity at low wind speeds. This could pos- 
sibly be corrected in the model by selecting a higher initial value for 
the drag coefficient. 

E. RESULTS FROM INDEX OF REFRACTION COMPUTATIONS 

A final application of the boundary layer model was the computa- 
tion of the index of refraction structure constant (0^) using (84) . The 
ratio of the exchange coefficients was determined from (85) and the 
input data for this equation were obtained by including the stability 
effect, the thermal wind effect, the correction for curvature, and using 
k=0.35. All parameters needed to compute were determined in the 
model . 

The results from these computations appear in Figure 36 . Values 

2 

of range from 0.0 in regions of near neutral stratification to a maxi- 
mum of 0.181 in the region of maximum wind speeds. Comparisons of 
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these values with those measured by Portman et al. (1968) indicate that 

they are comparable in regions of low wind speed with 0-9 less than 

a s 

2 

1 degree. Although no definite relation between and wind speed or 

2 

9 -9 appears in Figure 36, larger values of C,„ appear to be associ- 
a s T 

ated with the greater wind speeds and larger instability (9-0 ranging 

a s 

from -1 to -6°C). 

2 

The results indicate that a horizontal distribution of for a 
given height, in this case 19.5 meters, can be computed with the model 
with some degree of confidence in the range of values obtained. Com- 
parisons should be made between results computed by this model and 

2 

observations of to obtain better verification. 
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VI. CONCLUSIONS 



The results of this study have indicated that a marine boundary 
layer model can be applied to synoptic scale data to obtain various sur- 
face parameters. These parameters include friction velocity (U^), 
surface roughness (z^) , surface Rossby number (Rq), 19.5 meter wind 
(U ), veering angle (^^„)/ and stability length (z/L) , 

X b/ « O ^ 

First, the results indicate that the curvature correction produced 
the greatest change to the wind analysis. However, stability and modi- 
fication of the von Karman constant produced significant changes in the 
field, while changes due to the thermal wind were confined primarily to 
changes in wind direction. 

Second, the 19.5 meter wind was computed including various 
physical specifications and corrections and compared with synoptic 
observations. The results indicate that the computed wind was lower 
than the observed wind in all areas and for all cases, but the general 
patterns agreed well with observations . 

Finally, it appears that Card one's planetary boundary model pro- 
vides a reasonable framework to specify those surface layer parameters 
required for spectral wave forecasting and index of refraction calculations. 

Further studies are needed to more completely verify and improve 
the model. The specification of the drag coefficient is probably the 
aspect in this model which needs the most attention. The drag coefficient 
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has received considerable examination in observational investigations 
and the results indicate some disagreement among investigators. A 
suggestion of this study is to increase the initial value of to increase 
the computed surface wind speeds especially in areas of low wind speed. 
Application of this model to a region of reliable synoptic data such as 
during the BOMEX experiment in 1969 would reduce some of the errors 
due to data input and also test the physics of the model in the tropical 
atmosphere. Finally, observations of the index of refraction over a wide 
range of wind speed and stability would allow reliable comparisons to be 
made of the structure constant (C'j). 
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FT5ST = SIN(15.-"RADIAN) 

IF(S!NPHI.LT.FTSST) F = 2.*=0MHGA*FT£ST 
TH£RM = (GRAV*TOPT) /(F*>^-2‘MTL’MP (I ,J)+273. ) ) 
RHO = PR^;S(I,J )/ (RRMTFMPdtJ )+273. ) ) 

TD = OIF (I , J) 
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COMPUTE GEQSTROPHIC WIND FROM SFC P FIELD 
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ABSOLUTE VALU^: OF AIR-SEA TEMPERATURE DIFFERENCE LE 
CONSTANT DRAG COEFFICIENT = G.02A5 
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SUBROUTINE PSI(P,B,C) 
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